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Abstract 

We study magnetic fields in the charged condensate that we have previ- 
ously argued should be present in helium-core white dwarf stars. We show 
that below a certain critical value the magnetic field is entirely expelled from 
the condensate, while for larger values it penetrates the condensate within 
flux-tubes that are similar to Abrikosov vortex lines; yet higher fields lead to 
the disruption of the condensate. We find the solution for the vortex lines 
in both relativistic and nonrelativistic theories that exhibit the charged con- 
densation. We calculate the energy density of the vortex solution and the 
values of the critical magnetic fields. The minimum magnetic field required 
for vortices to penetrate the helium white dwarf cores ranges from roughly 
10 7 to 10 9 Gauss. Fields of this strength have been observed in white dwarfs. 
We also calculate the London magnetic field due to the rotation of a dwarf 
star and show that its value is rather small. 



1 Introduction and Summary 



Consider a system of electrically charged massive scalars and massive fermions of 
the opposite charge at high densities. When the interparticle separation is small, 
or the temperature of the system is high, neutral fermion-boson bound states will 
be unable to form. At high temperatures the equilibrium state of the system is a 
plasma. As the system cools below a certain critical temperature, the energy of the 
Coulomb interactions will significantly exceed the thermal energy, and the ionized 
system can crystallize. However, in certain cases the de Broglie wavelengths of the 
scalars begin to overlap before the crystallization temperature is reached. In this 
case, instead of crystallizing, the quantum-mechanical probabilistic "attraction" of 
the bosons can force the scalars to undergo condensation into a zero-momentum 
macroscopic state of large occupation number. The scalars minimize their kinetic 
energy, while phonons can not be thermally excited since the phonon mass gap 
- produced by this condensate - happens to be greater than the corresponding 
temperature. Therefore, after the phase transition all the thermal energy is stored 
in the near-the-Fermi-surface gapless excitations of quasi-fermions. We refer to this 
state as a charged condensate [1, 2, 3]. 

This condensation mechanism is different from that of the abelian Higgs model, 
or equivalently, the relativistic Ginzburg-Landau theory of superconductivity - the 
scalar field in our case has a conventional positive-sign mass term. A nonzero ex- 
pectation value for the electric potential (A ), or a nonzero chemical potential for 
the scalars, plays the role of the tachyonic mass, enabling the scalar field to acquire 
a vacuum expectation value [4]. The condensation mechanism is generic: the elec- 
tromagnetic interaction can easily be generalized to any U(l) abelian interaction, 
and the scalar field could be a fundamental field or a composite state. 

In this work we argue that the charged condensate has properties somewhat 
similar to type II superconductors. In particular, we show that it can admit solutions 
that are similar to the Abrikosov vortices [5], originally found in the Ginzburg- 
Landau model of superconductivity, and later recovered in the relativistic abelian 
Higgs model in [6]. The vortex solution is a topologically stable configuration, 
characterized by a nonzero winding number of the phase of the complex scalar 
field. Asymptotically, the scalar field is given by <p ~ ve ld , where v is the vacuum 
expectation value of the field and 9 is the azimuthal coordinate. Like the Abrikosov 
vortex, the charged condensate vortex carries a quantized magnetic flux. The vortex 
solution has a higher energy density than the pure condensate solution. However, in 
the presence of a sufficiently high external magnetic field, it becomes energetically 
favorable for the charged condensate to form vortices. 

The obtained vortex line solution exhibits the following structure: it has a narrow 
cylindrical core where the scalar field changes significantly from a nonzero to a 
zero value; this core is surrounded by a broad halo in which the magnetic flux is 
confined. The width of the latter region is determined by the penetration depth 
(i.e., the photon Compton wavelength). We refer to the system of the core and the 
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halo as the flux-tube, or the vortex line. This structure is similar to that of the 
Abrikosov solution. However, unlike the latter, our solution also carries a profile 
of the electrostatic potential within the halo, while this potential is exponentially 
small outside of the halo (i.e., the flux tube is charge neutral). Hence, in terms of 
differential equations, one has to solve three coupled equations instead of the two 
required in the Abrikosov [5], or Nielsen-Olesen cases [6]. 

One application of the charged condensate is to helium-core white dwarf stars. 
The cores of these dwarfs are composed of a highly dense system of helium-4 nuclei 
and electrons. At high temperatures, the electrons and nuclei form a plasma. While 
white dwarf stars composed of carbon, oxygen or heavier elements are expected to 
crystallize as they cool [7], it was argued in Refs. [2, 3] that the helium-4 nuclei would 
instead form the charged condensate, as the condensation temperature in the helium 
dwarfs is higher than the temperature at which the system would crystallize. This 
transition dramatically affects the cooling history of the helium-core white dwarfs. 
In fact, they cool faster; as a result, the luminosity function exhibits a sharp drop- 
off below the condensation temperature [3]. Such a termination in the luminosity 
function may have already been observed in a sequence of the 24 helium-core white 
dwarf candidates seen in NGC 6397 [8]. 

The above conclusions were obtained by considering white dwarfs that are not 
magnetized. Magnetized helium-core dwarf stars are also believed to exist; in mag- 
netized white dwarf stars, surface magnetic fields have been detected ranging from 
10 3 to 10 9 Gauss. If this is the case, it is important to know how the presence of a 
magnetic field would affect the above- described properties of the charged condensate. 

In analogy with type II superconductors, we would expect an external magnetic 
field to be entirely expelled from the charged condensate below a certain critical 
value of the field H c \. Above this value, however, we'd expect to have a mixed 
phase in which the magnetic field penetrates the charged condensate only in the 
form of the Abrikosov-like vortices. Finally, above a certain H c2 > H c i the magnetic 
field should entirely destroy the charged condensate. 

Indeed, this is the pattern that we find in the present work. The corresponding 
values of the critical magnetic fields in the interior of helium-core white dwarfs are: 
H c i ~ (10 7 — 10 9 ) Gauss, and H c2 ~ (10 13 — 10 15 ) Gauss, while in between these two 
scales we find quantized magnetic vortices permeating the bulk of helium-core dwarf 
stars. Hence in most of the magnetized helium white dwarfs the magnetic field will 
be expelled from the core in which the charged condensation has taken place. It's 
only the highly magnetized dwarfs, with fields ~ (10 7 — 10 9 ) Gauss, that can admit 
vortices. Also, since the value of H c2 ~ (10 13 — 10 15 ) Gauss is much higher that any 
magnetic field that may be present in dwarf stars, there will be no disruption of the 
charged condensate due to magnetic fields. 

The presence of the magnetic field would decrease the critical temperature at 
which the charged condensation would take place. However, this decrease will be 
significant only for the fields H close to the critical ones. For smaller fields, the 
change in the value of the critical temperature will be small AT C ~ —T C (H/H C ). 
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The organization of this paper is as follows: In section 2 we briefly describe the 
condensation mechanism for a generic system of scalars and fermions in the context 
of relativistic field theory. In section 3 we fix the phase of the scalar field to be 
of the vortex-type and solve the corresponding equations of motion. We compare 
our solutions to those found in the abelian Higgs model. In section 4 we consider 
corrections to our solutions due to the dynamics of the fermions. In section 5 we 
consider the effects of an external magnetic field on the charged condensate and 
determine the magnitude of the external field for which it becomes energetically 
favorable to form vortices. In section 6 we treat specifically the case of helium-4 
nuclei and electrons. We describe the system in the context of a non-relativistic 
low energy effective field theory, rather than the relativistic field theory used in 
the previous sections. We discuss the applicability of the vortex solutions found in 
section 4 to the helium-4 nuclei and electron system. We also consider the effect of 
a constant rotation on the condensate of helium-4 nuclei. 

2 Charged condensate: relativistic Lagrangian 

We start by considering a generic, highly dense system of charged, massive scalars 
and oppositely charged fermions. We assign a charge of +2e to the scalars and — e 
to the fermions in anticipation of the helium-4 nuclei and electron system to be 
discussed later in section 6. However, for now we keep our considerations general. 
The scalar field considered below could be any fundamental scalar field, possibly 
originating in beyond-the-standard-model physics. Our conclusions are independent 
of the specific charge assignment. 

The scalar condensate is described by the order parameter 0. A nonzero vacuum 
expectation value of implies that the scalars are in the condensate phase. Here we 
adopt a relativistic Lorentz-invariant Lagrangian which contains the charged scalar 
field 0, and the photon field A^. The fermion current is given by J M . 



The covariant derivative for the scalars is defined as = — 2ieA fl . 

The above Lagrangian could also contain a quartic interaction A(0*0) 2 . In the 
case that m% S> A J , this term will not alter significantly our results. Other possible 
terms including a Yukawa interaction between fermions and scalars were considered 
in Ref. [1]. For now we treat the fermions as a fixed, constant background charge 
density J M = Jo^o- We will relax this assumption in section 4 and consider effects 
due to the dynamics of the fermions and quantum loops of relativistic fermions. For 
simplicity we take the system to be at zero temperature (for some discussions of 
finite temperature effects see, [3]). 

Because the system has a conserved scalar current, we can associate with it a 
chemical potential fi s . For the Hamiltonian density, the inclusion of a chemical 
potential for the scalars results in the shift H — > H' = H — /i s Jg calar , where J^ calar = 




(1) 
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— i[(D &)*§ — <&*(Do<&)] is the time component of the conserved scalar current. For 
the Lagrangian density this shift can be written as a shift in the covariant derivative 
_D M — > D' — D^ — ifigd^o- In what follows primed variables H', CJ will refer to those 
variables which include a nonzero chemical potential for the scalars. 

The complex order parameter can be written in terms of a modulus and a 
phase = ^cre iCt . In these variables the Lagrangian density becomes 

v = -\K> + W^) 2 + K 2 ^ + - d * a ) 2 ° 2 - >h ° 2 - eA "J» ■ (2) 

From this form of the Lagrangian it is evident that a nonzero expectation value for 
Ao or a nonzero chemical potential fi s act as a tachyonic mass for the scalars [4]. 
In particular, when (2eA ) + fj, s — m H , the scalar field condenses, as we shall now 
show. 

Varying the Lagrangian with respect to and a gives the following equations 
of motion: 

—d^F^ = 2e{2eA v + fi s 5 u0 - d u a)a 2 - e J v , (3) 
Da = [(2eA li + ^ s 5^-d^af -m 2 H ]a . (4) 

Varying with respect to a gives the conservation of the scalar current: 

jscala, = + ^ _ a 2} = Q (g) 

On equation of motion (3) this expression is automatically satisfied. 

We now work in the unitary gauge where the phase of the scalar field is set to 
zero: a = 0. Note that this gauge choice is acceptable for a classical description of 
the condensate, however, in subsequent sections we will not be able to choose this 
gauge for the vortex solution for a well-known reason: in order to take a = 9 to 0, 
the corresponding gauge transformation A^ — > A^ + 8^9 would be singular at the 
origin where the vortex core is located. 

In the unitary gauge the equations of motion (3), (4) have the following static 
solution with a nonzero expectation value for a: 



(2eA )+fi s = m H , ( a ) = ^ . (6) 

In the condensate the gauge symmetry is spontaneously broken. The photon 
acquires a mass 



2e A /-^-. (7) 
V 2m H v ; 

This mechanism of symmetry breaking differs from the abelian Higgs model in that 
here the scalar field has a conventional positive-sign mass term. A nonzero expecta- 
tion value for A or a nonzero chemical potential fi s act as a tachyonic mass term. If 
we consider a system with no net charge then (A ) = 0. Then, from expression (6), 
in order for the scalar field to condense the chemical potential must satisfy ji s = ran- 
The bulk of the condensate is always electrically neutral, the scalar charge density 
exactly canceling the fermion charge density: 2eJQ Calar = ((2eAo) + /i s )(<r 2 ) = eJo. 



4 



3 Vortices in charged condensate 



To find the charged condensate solution (6), we fixed the phase of the scalar field to 
zero. We now consider a configuration where the phase is not set to zero, nor can it 
be set to zero everywhere by a non-singular gauge transformation. The requirement 
that the scalar field be single-valued everywhere is satisfied by demanding that the 
change in phase around a closed loop be an integer multiple of 2n. In a system with 
cylindrical symmetry, this is satisfied by setting a = n0, where 9 is the azimuthal 
coordinate and n is an integer. This phase can be removed everywhere by a gauge 
transformation — > A ^ + 8^(118), except at the origin where the scalar VEV goes to 
zero and the gauge transformation would be singular. The solutions of the equations 
of motion (3), (4) where the phase is fixed to a = n9 are vortex-type solutions. 

Hence, a = at the origin r = (where r is the 2D radial coordinate). Far from 
the origin however, we expect the solutions to recover the condensate solutions (6). 
At large r then, the gauge field takes the form 2eAj — > dja, or equivalently Ag — > 
n/(2er). From this form of the vector potential, it follows that this configuration 
has a quantized magnetic flux <P that is related to the integral of Aj around a closed 
loop at infinity: 



The magnetic flux is quantized in units of n. The quantization of flux implies the 
stability of the vortex configuration, although it may be possible for a high n vortex 
to decay into multiple vortices of smaller n. 

To solve the equations of motion (3), (4) for the vortex configuration we switch 
notation to dimensionless variables. The resulting equations are governed by a single 
parameter k, the ratio of the mass of the scalar to the mass of the photon in the 
condensate: k = m H /m 7 . This parameter k is the equivalent to the Ginzburg- 
Landau parameter in the theory of superconductors which gives the ratio of the 
penetration depth to the coherence length. For the helium white dwarf star, if we 
take the mass of the helium-4 nuclei to be roughly rnn = 3.7 GeV and the electron 
density to be Jo ~ (0.15 — 0.5 MeV) 3 , then we have n ~ 10 6 . In our derivations 
below we frequently take the large k limit. 

We define x = m 7 r, set A r = A z = and /i s = m H , and perform the following 
change of variables: 




(8) 



m^A^x) 
m y F(x) 
m H B(x) 



2exA g (x) , 

2ea(x) , 

/i s + 2eA (x) 



(9) 
(10) 

(11) 
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In terms of these new variables equations of motion (3), (4) become 



x — 



1 d 

1 d 

x dx 



dx 



d 




n 



x 



-M-lf F, 



2 



(14) 



(12) 



(13) 



The boundary conditions are set by requiring that the solutions asymptote to the 
condensate solutions for large r, while for r = we have A e = a = dA /dr = 0: 



We can compare these expressions to those obtained in the usual abelian Higgs 
model. Suppose that instead of Lagrangian (2) we had the abelian Higgs Lagrangian: 



£ AH = ~F% + \{d,af + \{2eA, - d,af a 2 v^f . (16) 



Using the same change of variables as above and defining m^ H = V\v : m^ K = 2ev, 
the equations of motion are: 



The equation of motion for the vector potential, expressed via A, is the same as in 
the charged condensate model. In the equation for the scalar field, the a 4 term in the 
abelian Higgs model gets replaced in the charged condensate model by a term that 
depends on the electric potential. In addition, in the charged condensate equations 
the electric potential is generally not zero and not constant and has its own equation 
to satisfy. 

Let us first examine the asymptotic behavior of the solutions to the condensate 
equations for x — > oo. Far from the origin we expect the fields to be very close 
to their condensate values. Then, on the r.h.s. of equation (12), it follows that 
A{x) — 1 = a(x) is very small. If we consider this equation only to first order in 
small fields then we can approximate the scalar field on the r.h.s. of (12) as F ~ 1. 
The solution for A that obeys the appropriate boundary conditions is 



For x -> : A(x) -> 0, F(x) -> 0, g -> . 
For x -> oo : A(x) -> 1, F(x) -> 1, B(x) -> 1 . 



(15) 




(17) 



(18) 



A(x) = 1 — c a xKi(x) . 



(19) 
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Here c a is a constant to be determined by the matching of the solutions and Ki(x) 
is the modified Bessel function of the second kind. In the large x limit this solution 
for A becomes 

A(x) ^ 1 - c^e'* . (20) 

To find the asymptotic behavior of B and F we expand these fields in terms of 
perturbations above the condensate values, B{x) = 1 + b(x) and F(x) = 1 + f(x), 
and we assume that b(x),f(x) <C 1. We then substitute these expressions as well 
as expression (20) into the equations for B and F and keep only terms linear in 
the perturbations b(x) and f(x). These two equations can be combined to obtain 
a fourth order differential equation for b(x). Using the ansatz b(x) = CbX s e~ kx 
where q>, s, and k are as yet undetermined constants, we can find the particular and 
homogeneous solutions for b(x) in the large x limit. We also take c a ~ 1 which we 
will justify later. For the particular solution we find 

, , \ Tin 2 e~ 2x 

U X) = M*TT) — - (21) 

For the homogeneous solution we deduce s = —1/2 and k 2 = (1 ± y/l — 1Qk 2 )/2. In 
the limit that k is very large, the solution becomes 

bh( x ) = — t=- \cism(y/~K,x) +c 2 cos(v / ^^)l , (22) 
with some constants C\ and c 2 . The complete solution is then 

B(x) = 1 + — - — — 1 [ci sm(y/K,x) + c 2 cos(y/Kx)] . (23) 

The solution for F(x) can be found once B(x) is known: 

= 1 + 8(/t 2 + 3) ~ 7^ ^ cos ^ a:; ) ~ 02 sin (V^ x )] ■ (24) 

Here c\ and c 2 are the same integration constants that appear in the expression for 
B(x). 

As k is large, the second term in the above expressions for B and F dominates 
the asymptotic behavior. For x^oowe have 

2 —It 
77-77 ^ p ^ x 

B W- 1 + 4RT3) — • (25) 

^)- 1 + 8^3)?- < 26 » 
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The asymptotic behavior for the vector potential and the scalar field are similar to 
that for the abelian Higgs model: 





c 



\/2nx 



(x) = l- C, 



a 



—x 



F (x) = 1 



(27) 



The vector potential A, and thus the magnetic field, are the same in both the charged 
condensate and abelian Higgs models. The asymptotic behavior of the scalar field 
in the abelian Higgs model in the large k limit is dominated by the e~ 2x term, as in 
the charged condensate model. 

Notably, this is not the asymptotic behavior for the Abrikosov vortex given in the 
Nielsen-Olsen paper [6]. This discrepancy was first pointed out by L. Perivolaropou- 
los in [9]. The incorrect asymptotic behavior is obtained if one similarly expands 
A(x) as A(x) = 1 + a(x) and only keeps terms linear in a(x). This is because the 
last term in (13) and the last term in (18) are quadratic in a(x) and yet, due to 
the different exponential dependence of the perturbations a(x), b(x) and f(x), these 
terms can be dominant over terms which are linear in b(x) and f(x). Linearizing A 
gives the correct asymptotic behavior of the fields only in the limit that k is small. 

The second term in the full expressions for B and F and thus the asymptotic 
behavior of both fields is due strictly to the presence of a nonzero magnetic field. 
In the absence of any magnetic field, the screening of any small perturbation of the 
fields above their condensate values vanishes as e~^ x = e~v" mHm ~< r (see Ref. [2]). 

The above asymptotic expressions are valid as long as x ^> 1, or equivalently 
r ^> 1/m.j. We consider now the intermediate region l/y/ft < i < 1, or equivalently 
1/M <C r -C l/m 7 where for convenience we have defined M = y/rrvfpm^. At 
distances much larger than 1/Mwe assume that the scalar field F is still close to 
its condensate value. Thus expression (19) is still valid for A. In this regime then 
n 2 /x 2 (A — l) 2 ~ n 2 /x 2 . The equations for B and F become 



Alternatively, we can once again expand B and F above their condensate values, 
B(x) = 1 + b(x) and F(x) = 1 + f(x), and solve for b(x) and f(x). The homoge- 
neous solutions for b(x) and f(x) are the same as those given above with the same 
coefficients C\ and c 2 . Solving for the particular solutions gives the full solutions in 




(28) 



The solutions are straightforward to find: 




(29) 
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the linearized approximation: 

B(x) = 1 + 2 H [ci sin(v/fi;a;) + c 2 cos(\/Ka;)] , (30) 

= 1 " 4^ + ^4 - —^f- h cos(^x) - c 2 sin(v^)] . (31) 

The coefficients C\ and c 2 are needed to perform the matching. However, as we'll 
see below, these coefficients will turn out not to be exponentially large, and hence 
solutions (30) and (31) approximate well the solutions in (29). 

The approximations made to find both the homogeneous and particular solutions 
break down as x approaches 1/ ' \[k. Moreover, f(x) becomes of order 1 at x ~ l/\/ft 
and thus the linear approximation in general no longer holds below this scale. 

Finally, we'd like to solve in the r — > limit. Before we do so, however, we 
emphasize that validity of this procedure needs some justification. The interparticle 
separation is given by d oc J ^ 3 . This corresponds to x oc l//t 1//3 . At distances 
shorter than this x we expect that an effective field theory would break down and 
thus it would make little physical sense to solve the equations (12), (13) and (14) 
in this regime. Moreover, the scale 1/M is typically shorter than the interparticle 
separation d, hence, particles at these scales cannot in general be modeled by a 
smooth distribution. 

However, both fermions and bosons are in a condensate state in which the lo- 
cation of individual particles has uncertainties much greater than the interparticle 
separation. Hence the latter notion loses its meaning as a microscopic characteris- 
tic of the system. For this, we'll still approximate particle distributions by smooth 
functions all the way down to the scale ~ 1/M, which is a dynamically determined 
short-distance scale at which weakly coupled expansion breaks down [1]. As to solv- 
ing at scales less that 1/M, we regard this as a purely mathematical exercise aimed 
at finding the matching of the asymptotic solutions for the corresponding differential 
equations for all values of the coordinate x. 

Taking A, B, and F to be series expansions in small x obeying the appropriate 
boundary conditions, the solutions to (12), (13) and (14) are 

A(x) = a x 2 -^-x\ (32) 

o 

B(x) = 60-^, (33) 



F(x) = f 



x-\ {n 2 {bl - 1) + 2a ) x 3 

o 



(34) 



For simplicity we have solved for the case that the winding number n — l. 1 The 
coefficients a , fo, and b , as well as the coefficients c a , c±, and c 2 can be determined 



1 For ti/ 1, the leading term in the expansion for F will be oc x n . The expression for B{x) 
remains unchanged and the leading term in the expansion for A(x) = a^x 2 is also the unchanged. 
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by matching the above solutions to those in the intermediate region, given by (19), 
(30) and (31). 

To determine the physically appropriate matching radius, we first use Gauss's 
law to find the charge of the vortex solution. The number density of fermions in the 
vortex Jo is constant and is the same as the number density of fermions in the normal 
condensate phase. We have fixed it so by hand, but will justify this later. The scalar 
number density is given by ^J BF 2 and varies as a function of x. Therefore it is 
not in general equal to its condensate value | J . The variation of the scalar number 
density away from its condensate value can lead to a net charge density of the vortex 
core within the vortex halo. In particular, there are two competing effects. In the 
intermediate region 1/v^ < i < 1, both B and F are above their condensate 
values, thus the scalar number density is greater than the scalar number density in 
the condensate. As x — > 0, however, F — > and the scalar number density drops 
to zero, significantly below the condensate value. The matching radius should be 
chosen so that these two effects combine to give the appropriate charge density as 
determined by Gauss's law. 

From Gauss's law we can calculate the the average charge density of the vortex 
inside radius x — 1. As is usually the case, we can determine the net charge enclosed 
in a region knowing only the form of the potential at the boundary of that region. 
Equation (29) gives the potential in the intermediate region independent of matching 
coefficients c\ and C2- This form of the potential, together with Gauss's law, allows 
us to calculate the net charge of the vortex at x — 1 independent of the matching 
conditions and the x — > solutions. 

Gauss's law is given by equation (3): 



The r.h.s. of the above equation is the charge density. Integrating both sides of the 
above expression over the volume of the vortex and dividing by the total volume 
gives the average charge density inside distance x: 



Here V is the volume equal to the length of the vortex times the cross-sectional area 
and the r.h.s. is evaluated at the boundary of the vortex. Using expression (29) 
for B(x) at x — 1, the average charge density inside x — 1 is Q enc /V = — 2eJ //t 2 . 
The negative sign indicates a dearth of scalars in this region, but, as k is very large, 
this is a small correction to the overall average charge density of scalars oc eJ . 
To check this result one can likewise use the asymptotic solution for B, expression 
(23), at x — 1. Assuming that the coefficients C\ and c 2 are not exponentially 
large and thus these terms are not dominant in the solution for B at x = 1, one 
finds Qcnc/V oc — 2eJo//t 2 . This is consistent with the previous result. Farther out, 
B(x) — 1 is exponentially suppressed thus the net charge of the vortex approaches 
zero as x becomes large. 




(35) 



Q. 



cue 




(36) 



V 
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We can now use this result to determine the matching radius R. Given the 
smallness of the average charge density found inside x — 1, the excess of scalars 
in the intermediate region of the vortex must cancel the shortage of scalars in the 
x — > region to great accuracy. Using expressions (29) for B and F, it can be shown 
that this happens when R ~ l/y/R. Thus we use this as our matching radius R in 
what follows. 

We start by matching the solution for A(x) in (32) and its first derivative with 
its solution in the intermediate region (19). Taking the matching radius to be small, 
gives: 



ao = ~2 



■2 



c a = l + ^, (37) 



where 7 is the Euler-Mascheroni constant. As long as R is less than one, ao is 
positive. Moreover, we see that we were justified in taking c a ~ 1 in our previous 
calculations. The magnetic field is given by 

ff=^i" (38) 
2e x ax 

Near the origin the magnetic field is of order m^/(2e). For x > 1/a/k it is given 
by rn^Ks(x) / (2e) , where Kq{x) is the modified Bessel function. For x 1, i.e. for 
r ^> l/m 7 , the magnetic field is exponentially small. 

To find the remaining coefficients, we use a obtained above and match (33) 
and (34) and their first derivatives to the appropriate solutions in the intermediate 
region (30), (31). We now take the matching radius to be R = The solution 

with the lowest energy is one in which the scalar field F(x) is identically zero in the 
region x < R. The corresponding coefficients are 

b = l + l, d = K - 5 / 4 e(2cos(l) + sin(l)) , 
f = 0, c 2 = /t- 5 / 4 e(cos(l) -2sin(l)) . 

In Fig.l below, the fields are plotted for small r and for k ~ 10 6 . The radius r = 1/M 
corresponds to the matching radius x — R — 1/ yJH. The radius r = d denotes the 

— 1/3 

interparticle separation d — J . Unlike the magnetic field, the potential and 
scalar field approach their condensate values for x > l/^fn. This is in contrast to 
the abelian Higgs model in which the scalar field is close to its condensate value for 
x > 1/k, i.e. for r > l/mn- 



4 Fermion Dynamics 

In our discussions above we have treated the fermions as a fixed charge background 
J pl = Jq^q. We relax this assumption now and introduce dynamics for the fermions 
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F(r) B(r)-] 




(a) Scalar field as a function of radius (b) Potential as a function or radius 

Figure 1: Small r solutions for the scalar field and electric potential 

via the Thomas-Fermi approximation. The fermion dynamics are governed by the 
constant chemical potential hf'- 



\i>F = \JPf( x ) +™ 2 e + eA (x) . (40) 

The local number density of fermions is determined by the Fermi momentum: 
Jq(x) = p F (x) I (Stc 2 ) . In this way the number density of the fermions Jo gets related 
to the electric potential A . For relativistic fermions 

Jo(x) ^ T^(fi F - eA (x)) 3 . (41) 

The chemical potential gets fixed by the value of the fermion number density in the 
condensate phase, where (Aq) = 0. If Jo represents the number density of fermions 
in the condensate, then fjip — (37r 2 J ) 1//3 . The photon mass m 7 is also defined in 
terms of J : m 7 = 2e\J J /2rriH- In the vortex phase Jq{x) — > Jo for large x. 

To include the effects of an x-dependent Jo into our equations, (41) gets incorpo- 
rated into the equations of motion (3). As a result the equation of motion for B(x) 
(14) gets modified. In the linearized equations, the effect is the addition of a new 
term for b(x) with a coefficient which scales as cx tjih/^f- However, it turns out that 
this new term does not contribute significantly to the solutions. This is because, in 
the fourth order differential equation for b(x), terms with coefficient itih/^f oc k 2//3 
are subdominant compared to terms with coefficient k 2 . Accordingly, the solutions 
found above in the intermediate and asymptotically large regions are still valid. It 
can be shown that the x — > solutions (32), (33), (34) are also unaffected. In 
physical terms, the inclusion of the fermion dynamics via the Thomas-Fermi ap- 
proximation gives rise to ordinary Debye screening. This screening is subdominant 
compared to other screening effects in the condensate (see [2]). Moreover, the profile 
of A (x) away from the core and within the halo is very shallow, giving rise to a very 
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mild dependence of the charge density on x. Hence, the latter can be approximated 
by a constant, as was done in the previous sections. 

The Thomas-Fermi approximation does not capture the possibility of exciting 
gap-less modes near the Fermi surface. To include this effect we must calculate 
the one-loop correction to the gauge boson propagator. In other words, we must 
include in the Lagrangian (1) the fermion kinetic, mass and chemical potential terms 
and take into account the known one-loop gauge boson polarization diagram when 
calculating the gauge boson propagator. This was done in the second reference in 
[2] and we use those results here in what follows. 

The one-loop correction to the gauge boson propagator gives corrections to the 
static potential A . We are interested in how this correction compares to the poten- 
tial found in the intermediate region of the vortex (30). To estimate its magnitude, 
we consider a toy model of the vortex. We find the potential due to a wire of con- 
stant linear charge density A located at r = 0. The linear charge density of this 
wire is set by the characteristic charge of the vortex: since the scalar charge density 
varies significantly from its condensate value eJo at scales r < 1/M, it follows that 
at short distances the linear charge density of the vortex can be approximated by 

At large distances the vortex is effectively neutral, as mention above. Thus we expect 
the one-loop contribution to the static potential to be irrelevant at large scales. 
In three dimensions, the charge density of the source is given by 

J sourcc (r, 6, z) = —5(r) . (43) 
7rr 

The static potential is determined from this source and from the {00} component 
of the gauge boson propagator D 00 : 



A (r) = - J d 3 r' D 00 (r - r') J sourcc (r') . (44) 
The propagator was found in the second reference in [2]: 



-£>oo(w = 0,k) = ( k 2 + m 2 + — + F(k 2 ,k F ,m f )) 



-i 



(45) 



The function F(k 2 , k F ,rrif) is due to the one-loop photon polarization diagram. 
It includes both the vacuum and fermion matter contributions. Here kp denotes 
the Fermi momentum and rrif the mass of the fermion. A complete expression for 
F(k 2 ,k F ,rrif) can be found in Ref. [10]. We take the expression for F(k 2 ,k F ,rrif) 
in the massless (rrif — 0) limit, which is a good approximation for ultra-relativistic 
fermions: 
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Here /xq stands for the normalization point that appears in the one-loop vacuum 

polarization diagram calculation. 

The function F introduces a shift of the pole in the propagator, corresponding 

to the "electric mass" of the photon. This part of the pole can be incorporated 

via the Thomas- Fermi approximation, as was done above. In addition, however, the 

function F also gives rise to branch cuts in the complex |k| plane (see Ref.fll] for the 

list of earlier references on this). These branch cuts give rise to the additional terms 

in the static potential which are not exponentially suppressed, but instead have an 

oscillatory behavior with a power-like decaying envelope. In a non-relativistic theory 

they're known as the Friedel oscillations [11]. In the relativistic theory they were 

calculated in Refs. [12, 10]. We have calculated them in Ref. [2] for the relativistic 

theory in the presence of the condensate: taking the Fourier transform of (45), the 

dominant contribution due to the branch cuts is 

_ a cm k 5 F sm{2k F f) 
-D m {r) - M8 _ 4 . (47) 

Here f represents the 3D radius in spherical coordinates, as opposed to the 2D radius 
r. 

Using this expression together with expression (43) in equation (44), the correc- 
tion to the static potential is 

a cm X k F f°° , sin(2 + r 2 ) 
Mr) = ^^pL dZ (z' 2 + r 2 ) 2 • (48) 



An upper bound on the potential can be found by taking sm{2k F \/ z' 2 + r 2 ) 
After integrating, this gives 



A (r) < -^-^ttti TV 2~T ( 49 ) 

7i M 2r 3 e 2 V m H m 2 H r 3 

On the vortex solution, the leading term in the potential in the intermediate 
region is given by expression (30): 

Mr) = ^(B{r) - 1) ~ —^— 2 . (50) 
2e 4em#r z 

Given that both kplvriu 1 and \j{mnr) < 1, we see that the one-loop correction 
to the potential (49) is greatly suppressed compared to the potential found in the 
vortex solution. Thus the excitations of the fermions do not significantly alter the 
vortex solutions. 

One further effect that we take into consideration is the Landau quantization of 
the fermion energy levels due to the presence of the magnetic field in the interior 
of the vortex. In the presence of an applied magnetic field, the separation between 
the Landau levels is given by uj = eH/rrif where H is the magnetic field and m/ is 
the fermion mass. Near the core of the vortex where H ~ m 2 /(2e) the separation of 
levels of the fermions is uj ~ m 2 /mf. Since the photon mass m 7 is generally much 
smaller than the fermion mass, this shift in energy is negligible compared to the 
typical energy of the fermions. 
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5 Energetics and external fields 



We now turn to the question of when it is energetically favorable to form a vortex 
in the charged condensate. We start by comparing the average energy density of 
the vortex to the energy density of the pure condensate. Above, using Gauss's law, 
we found that inside the distance x — 1 the vortex has a small negative charge 
density implying that in this region the average scalar number density is lower than 
in the condensate phase. At distances x ^> 1 this charge density is exponentially 
suppressed indicating that the net charge of the vortex is zero and thus the total av- 
erage scalar number density is the same in both the vortex phase and the condensate 
phase. In calculating the average energy density of the vortex inside the distance 
x — 1, we are not interested in the contribution to the energy due to the discrep- 
ancy in the number of scalars between the vortex phase and the condensate. This 
contribution to the overall difference in energy vanishes at large distances. Thus we 
calculate the energy density of the system using TV — Ji — fi s jQ Calar . The additional 
term effectively subtracts off the energy density due to the scalar number density. 
We compare H' in the vortex phase to W in the condensate. 

The Hamiltonian density H' can be calculated from the Lagrangian CJ (2): 



We have simplified the Hamiltonian using the equations of motion (3), (4) and have 
taken boundary terms to be negligible. The magnetic field H and the electric field 
E are defined as usual 



The fourth term in the Hamiltonian is exactly — h s Jq r as we would expect. The 
third term is due to the energy of the scalar field. Unlike in the abelian Higgs model, 
the energy density of the scalar field in the center of the vortex, i.e. in the "normal" 
phase, is lower than in the condensate phase. However, the energy density of the 
vortex is still greater than that of the condensate alone, due to the gradients of 
the scalar field and due to the intermediate region l/y/n < i C 1 in which the 
values of both the potential and the scalar field are greater than their condensate 
values. This contribution to the energy density is roughly equal in magnitude to 
the contributions coming from the electric and magnetic fields. On the condensate 
solution, the Hamiltonian density is identically zero: 7~C' CC = 0. 

For large x, deviations away from the condensate are exponentially suppressed 
and thus differences in energy between the two phases are negligible. So to find 
the average energy density within the vortex, we integrate the Hamiltonian density 
over an area of radius x = 1 and then divide by the total area. The average energy 
density within the radius x < 1 is 



W = \H 2 + \E 2 + |(2eA) + n 9 - a) 2 a 2 - fi s (2eA + /i s - a)a 2 



(51) 



H = ~(rA (f ), E = 
r dr 



dr 



-o 



(52) 




(53) 
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The Hamiltonian density can be further simplified using the equations of motion 
(see Ref. [13] for more details). In terms of the dimensionless variables defined 
above equation (53) becomes: 

i r 1 

e ave = ^n H J J xdx (B(x) - 1) (3 - F{xf) . (54) 

For the region x < 1/ \/k we use solutions (33) and (34) and in the intermediate 
region 1/y/it < x < 1 we use solutions (30) and (31) with the coefficients found from 
matching (39). Upon integration, the average energy density is 

eave = ^(logK+14) . (55) 

The numerical coefficients should not be taken too literally given the approximations 
made in obtaining the solutions which yield the above result. However, the overall 
scaling of the energy density e ave oc m H J {\ogK) / ' k 2 is remarkably independent of 
the matching radius and other details of the solutions. As the energy density of 
the condensate is effectively zero (J~i'cc = 0); the aDove expression represents the 
difference in energy between the two phases. 

To see when it is energetically favorable for the condensate to form vortices, we 
now consider placing the condensate in an external field H ext pointed along the z- 
axis. We shall see that the magnetic properties of the charged condensate resemble 
those of a superconductor. In particular, when k 1, the charged condensate 
resembles a type II superconductor. When an external magnetic field H cxt is applied 
to the condensate, below a critical value H c \ surface currents oppose the penetration 
of the field and the induction B ind is zero in the bulk of the condensate. For H cxt > 
H c i magnetic flux penetrates the condensate in the form of vortices. At another 
critical value of the magnetic field H C 2 the normal phase is restored and the induction 
B im \ is equal to the applied field H ext . In what follows we determine the critical values 
of the fields H c \ and H c2 . 

Given the energy density e of the vortex phase above the pure condensate phase, 
we can find the value of the magnetic field H cl at which it becomes energetically 
favorable to form vortices. In the absence of an external field, it is never energetically 
favorable to form vortices as the energy density of a vortex is greater than that of the 
pure condensate. In the presence of a small external magnetic field, below H c i, the 
condensate must expel the magnetic field entirely from its bulk in order to remain 
in the condensate phase. This requires energy; the energy per volume needed to 
expel the external field is ^H 2 xt . If vortices form in the condensate then the energy 
required to expel the magnetic field is smaller than if the field were to be completely 
expelled. More specifically, if the vortices give rise to an average magnetic field in the 
condensate B ind , then the energy needed to expel the remaining magnetic field would 
be |(i?ext — -Bind) 2 - The energy gained by forming vortices is the difference between 
this energy and the energy required to expel the magnetic field entirely. Assuming 
that And is small compared to H ext near the transition point, this difference can be 



16 



approximated by ^(2Bi n ^H ext ) . Thus for a high enough external field, the energy e 
lost in creating a vortex is compensated by the energy gained in expelling a smaller 
magnetic field B ind H cxt . In order for formation to be energetically possible, we must 
have e < B ind H ext . The equality determines the critical external field H c \. (See Ref. 

[5]-) 

Suppose the number of vortices per area in the condensate is given by N. Then 
the energy density due to the formation of vortices is given by e = iVA where A is 
the energy density per unit length of a vortex. Using e avc found above (54) as the 
energy density of a single vortex, 

7T 

A = -^-e avc . (56) 

The induction £> in d is given by 

B ind = N j A-dl=^. (57) 
Combining these expressions, the critical field H ci = N\/(Bi nd ) is given by 

o 

H cl = -1 (log(«) + 14) . (58) 
8e 

The final expression for H c \ is independent of the number density of vortices N. 
It follows that if it is energetically favorable to create one vortex, then it will be 
even more energetically favorable to create many, up to the point than interactions 
between vortices become significant. At distances greater than r = l/m 7 we expect 
fields outside the vortices to be exponentially suppressed and thus the vortices to 
be effectively non-interacting. So at the transition point H c i, it is likely that the 
number density of vortices is of the order iV ~ m 2 j 



7T. 



If we take Jo — (0.15 — 0.5 MeV) 3 , a reasonable value for white dwarfs, this 
gives a magnetic field of roughly H cl ~ (10 7 — 10 9 ) Gauss. Thus, the vortex lines 
should be expected to be present in the bulk of the helium-core white dwarf stars 
with strong enough magnetic fields. 

A sufficiently high magnetic field will disrupt the condensate entirely. One way 
to approximate the magnetic field at which this transition occurs is to consider the 
density of vortices in high external magnetic fields. When the cores of the vortices 
begin to overlap, then the scalars are mostly returned to the normal phase. We define 
the core of the vortex to correspond to x — 1/y/K, or equivalently, r = 1/M as this is 
the region in which the VEV of the scalar field drops to zero. As N is the number of 
vortices per area, at the transition point iV ~ M 2 /tt. We define the critical external 
field at this point to be H c2 . When the condensate enters the normal phase, the 
induction B ind is equal to the external magnetic field H c2 = B ind = 2nN/ (2e). From 
these two expressions we find H c2 : 

M 2 

H c2 = —. (59) 

e 
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For J ~ (0.15-0.5 MeV) 3 , H c2 ~ (10 13 -10 15 ) Gauss. This is well above the values 
of the fields expected to be present in a majority of white dwarf stars. Thus the 
external magnetic field is unlikely to be large enough to push the condensate into 
the normal phase. 

We note that both H c \ and H C 2 given above were determined at zero temperature. 
Generally, we expect these expressions (58), (59) to be valid at temperatures well 
below the condensation temperature. 

Finally, in type II superconductors the dependence of the critical temperature on 
the magnetic field is well-approximated by T c /2 /T c 2 ~ (H c — H)/H c , where T' c is the 
transition temperature when the magnetic field H is present. We expect a similar 
relation to be valid in our case too. Hence, as long as the value of the magnetic field 
is not too close to either critical value, the change of the transition temperature due 
to the magnetic field should be small. Near the critical values, however, the change 
of the phase transition temperatures (from the normal to the vortex phase and from 
the vortex phase to the phase with no magnetic field) could change significantly. The 
would be crystallization temperature will also change, and the charged condensation 
may or may not be favorable for close-to-critical magnetic fields 2 . 

6 The low energy effective Lagrangian 

For a system of helium-4 nuclei and electrons, one can consider a non-relativistic 
effective Lagrangian for the order parameter $, as the helium nuclei are non- 
relativistic in the condensate phase. This Lagrangian is less restrictive than the 
relativistic Lagrangian in that it is not required to be Lorentz invariant. The non- 
relativistic Lagrangian must give rise to the Schrodinger equation for the order pa- 
rameter in lowest order in the fields and it must respect the appropriate symmetries 
of the physical system, including translational, rotational and Galilean symmetries 
as well as gauge invariance. (See [2, 3] for more details.) Such a Lagrangian was ini- 
tially proposed by Greiter, Wilczek and Witten in the context of superconductivity 
[15]. We used it here to describe the charged condensate: 

Ceff = V (^(<t>*A)<f> - (A,*)**) - , (60) 

where D = (do — 2ieA ), Dj = (dj — 2ieAj), and V(x) is a polynomial function of its 
argument. We could introduce a chemical potential for the scalars into the argument 
of V(x) in the form fi^fi^*^. The relationship between the relativistic chemical 
potential fx s and the non-relativistic chemical potential is given by ji s = inn + (J>nr- 
Thus the neutral condensate where (A ) = and fi s = ran corresponds to unr = 
(assuming that the quartic and other interactions are neglected). Again, one could 

2 For discussions of other magnetic field effects in very highly magnetized white dwarfs, see [14] 
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also include a quartic term A($*$) 2 . However, as long as w? H ^> AJ and A < 1, 
this term can be neglected 3 . 

In the condensate where the VEV of $ is nonzero, we can express <£> in term of a 
modulus and phase: $ = £ exp(ir). Written in terms of fields £ and T, the effective 
Lagrangian (60) takes the following form: 

C eff = V (Wo - d m 2 - ^(V,S) 2 - 2^(2e^ - d,T) 2 ^ . (61) 

The gauge field couples to the electron density as — eA J , thus we include this term 
in the Lagrangian (61). Once again we work in the unitary gauge and set T = 0. The 
equations of motion which follow from (61) then have the following static solution: 

2e£ 2 = eJ , = 0, V'(0) = 1 . (62) 

The above solution describes a neutral system in which the helium-4 charge density 
2e£ 2 exactly cancels the electron charge density —eJ . Since on the solution the 
argument of (61) is zero, the condition V'(0) = 1 is satisfied by any polynomial 
function for which the first coefficient is normalized to one: V{x) = x + C 2 :r 2 + .... 

The condensate solution sets a preferred Lorentz frame. We consider the dy- 
namics of small perturbations above this background. We express £ in terms of a 
perturbation r above the condensate value: 



£(x)=vW| \l^ + r(x) I . (63) 



2m 



H 



The Lagrangian, including the gauge field kinetic term, expanded to second order 
in fields becomes: 



leff = - \id 3 rf + \c 2 m H J Q m 2 Al - \m 2 A) + 1m H m^r . (64) 



We can compare this Lagrangian to the one obtained from the relativistic theory 
(2) by enforcing Lorentz invariance, i.e. by demanding that the "electric" mass 
of the gauge field be equal to the "magnetic" mass. This would fix the value of 
C 2 = l/(m H J ). With this value of C 2 , the Lagrangian for small perturbations 
above the condensate in the non-relativistic theory is identical to the Lagrangian 
for small perturbations in the relativistic theory found in Ref. [1], up to a time 
derivative for r. 

However, we do not in general expect that the low energy effective theory will 
obey the Lorentz invariant condition C 2 = l/(m#Jo). Instead, C 2 must be fixed by 
the particular physics of the system. It's worth noting that even if C 2 were to be 
fixed via Lorentz invariance, introducing fermion dynamics via the Thomas-Fermi 



3 The chemical potential and quartic terms would have to be retained if we were to discuss 
temperatures near the phase transition point. 
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approximation introduces an additional term into the Lagrangian (64) of the form 
e 2 /j 2 J 3 Al/-K 2 which breaks the degeneracy between the electric and magnetic masses. 
This is a typical scale by which we'd expect the electric and magnetic masses squares 
to differ from each other. 

In section 3, in order to find vortex solutions in the intermediate region 1 / \/k <C 
i < 1 and the asymptotic region x ^> 1, we treated the electric potential A and the 
scalar field a in the linear approximation, but kept higher order terms for the vector 
potential. Thus to determine the applicability of the solutions found above to the 
non-relativistic effective theory, we should consider higher order terms in Aj than 
the ones given in (64). We also restore the phase T. Given V(x) = x + C 2 x 2 + 
the equations of motion to next-to-leading-order are 

-d^F^ = 2e [1 + 2C 2 S 2 (2eA - d T)} £ 2 - eJ , (65) 
—d^F^j = 2e(2eA j - o^E 2 , (66) 
-V 2 S = [2m H (2eA - d T) + AC 2 m H Y?{2eA Q - d T) 2 ] E - (2eA j ) 2 S . (67) 

If we take £ = ^/itlhU and C 2 = l/(m#Jo), then the first two equations of motion 
above (65), (66) are the same as in the non-relativistic case (3), up to second order 
in small fields. The third equation (67) has an extra factor of (2eA ) 2 S compared 
to equation (4). However, since this term is second order in A and we treated A 
in the linear approximation, this does not alter our solutions in the intermediate 
and asymptotic regions. Thus for C 2 = l/(m#Jo), the vortex solutions found above 
for x ^> 1/-\/k are also solutions for the non-relativistic effective theory. As we 
mentioned above, however, we would expect the realistic value of C 2 to be different 
from the one we used by the quantity (e 2 /j 2 J 3 / m H J m 2 ) . However, as we have shown 
in section 4, such corrections are subdominant because the value of the photon 
electric mass is smaller that the value set by the scale M. 

As in the relativistic case, the solutions formally break down near x = when 
the change in the scalar field becomes of order 1 and thus the linear approximation is 
no longer valid. More realistically though, we do not expect the effective field theory 
to hold at distances shorter than the interparticle separation x oc l//t 1//3 . Instead it 
will cease to be a valid description of the physics before reaching x — 1/ y/R. 

We can use the non-relativistic formalism to consider the effects of the rotation 
of a white dwarf star on the magnetic field in its interior 4 . In this formalism the 
scalar number density and current density are given respectively by: 

jscalar = ^ jscalar = [(£1^)*$ _ . (68) 

The number density is related to the current density by Jj calar = J^ caUr Vj where Vj 
is the velocity vector of the rotating scalar particles. Using the change of variables 



J We thank Daniel Stein for raising this issue. 
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defined above, $ = Eexp(ir), we can use these expressions to find vf 

Vj = ^ (2e^- - 9jF) . (69) 

This known result is notably different from that of a superfluid in which the scalar 
field does not couple to a gauge field. In the absence of the Aj term in the above 
expression, one would conclude that V x v = and thus the scalar condensate does 
not support rotation. Instead, in the presence of the gauge field we find 

V x v = — H . (70) 
m H 

The magnetic field H is called the London field [16] . 

The velocity vector v can be written in term of the angular velocity v = fl x r. It 
follows that, for constant fl, the rotation of v is given by V x v = 2fl. Accordingly, 
the magnetic field can be expressed in terms of the angular velocity: 

2e m 2 K J 

Here Jo is the fermionic number density. Thus the condensate of helium-4 nuclei 
can rotate with the rest of the star, unlike a neutral condensate. The consequence 
is a small, constant magnetic field in the bulk of the condensate. 

Varying the Lagrangian (61) with respect to Aj gives 2eJj calar = eJj, where Jj 
is the fermion current density. Using Jj = J$Vj, it follows that the fermion velocity 
vector is equal to the scalar velocity vector. The electrons and the helium-4 nuclei 
rotate together in the core of the star. At the surface however, there is a thin layer of 
helium-4 nuclei that is slightly out of rotation with the rest of the star. This feature 
becomes evident upon finite volume regularization of the system. The thickness of 
the layer is roughly l/m 7 . This surface layer is what gives rise to the London field 
in the interior of the star [17]. To estimate the value of the London field we take 
the angular velocity of a helium white dwarf star to be Q ~ 10~ 2 Hz. The resulting 
London field is H ~ 10~ 6 Gauss. This field is present even in the absence of vortices. 
However, it is too small to affect any of the results given above. 
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